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ОПТИМІЗАЦІЯ МАТЕМАТИЧНОГО І ПРОГРАМНОГО ЗАБЕЗПЕЧЕННЯ ДЛЯ 
ЧИСЕЛЬНОГО РОЗВ'ЯЗАННЯ СЛАР В ОБЕРНЕНИХ ЗАДАЧАХ ГРАВІМЕТРІЇ 

В розвитку сучасних засобів інтерпретації даних геофізики спостерігається стійка тенденція до створення спеціалі- 
зованих програмних пакетів та інтегрованих середовищ для швидкого й надійного комп'ютерного експрес-аналізу та тлума- 
чення цих даних майже одразу після вимірювань в польових умовах. В умовах зростаючих запитів щодо точності, швидкодії 
і зручності користування таких пакетів з'явились відомі програмні продукти як в ближньому (Сигма-30), Вектор, Согато- 
дог), так і дальньому (5игСЕ, СМ-5увз) зарубіжжі. Другий напрям в інтерпретації й моделюванні геофізичних даних передба- 
чає використання окремих програмних модулів з розвинутих бібліотек програмного забезпечення для задач обчислювальної 
математики (ІМ5Ї,, Н51., АТІ.А5, ВІ А5, МАС, Нагмеї!, тощо) з метою розв'язування систем лінійних алгебраїчних рівнянь 
(СЛАР), до яких зводяться в математичній постановці окремі математичні моделі геологічного середовища. Про успішне 
використання такого підходу на основі підпрограм з бібліотеки ІМЧІ, вказує робота |1). 

Крім того для моделювання, візуалізації й перевірки гіпотез досить широко використовують універсальні матема- 
тичні пакети (Маріє, Мафетайса, Мафар, МафСад, РДегіуе, МиРад, 5сії аб) та графічні інформаційні системи (ПаСео, 
АгсУіем/). Згадані вище пакети й бібліотеки пропонують розрізнений набір методів, як правило, обмежених за точністю або 
порядком вирішуваних систем (21, або ж "заточених" під геологічну будову й фізичні обмеження конкретних територій. Цих 
методів достатньо для вирішення загальних задач лінійної алгебри, але не достатньо для розв'язання специфічних задач, що 
виникають в обернених задачах граві- та магнітометрії, що вимагають регуляризації і багатокрокових стратегій ітераційного 
покращення розв'язку. До того ж, практично всі ці пакети комерційні 1 досить дорогі. З огляду на це актуальність розробки 
подібного вітчизняного програмного забезпечення не викликає сумніву. 

Задля стійкого чисельного розв'язання реальних задач структурної гравіметрії потрібно розвивати, як окремі модулі 
й надстройки, програмні засоби для вирішення СЛАР з розрідженими й погано обумовленими матрицями великої розмір- 
ності, застосовувати спеціальні методи апроксимації обернених операторів й оптимізації розв'язків. В згаданих вище засобах 
таких модулів немає. 

Нами планується розробити комплекс програм, адаптованих для сучасних комп'ютерів, який передбачає розв'язання 

систем лінійних алгебраїчних рівнянь з симетричними й несиметричними розрідженими матрицями, до яких зводяться обер- 
нені задачі граві- 1 магнітометрії. В цьому комплексі передбачається розробка наступних блоків: 
вирішення прямої задачі гравіметрії для контактної поверхні; 
вирішення оберненої контактної задачі для поля, заданого на короткому інтервалі; 
аналітичне продовження гравіполя в верхній і нижній півпростір; 
згладження й апроксимація даних виміряного гравіполя; 
перетворення вхідних даних в різних картографічних проекціях; 
вирішення систем лінійних рівнянь з розрідженими матрицями. 
Пряму задачу гравіметрії для кількох контактних поверхонь вирішуємо методом підбору шляхом апроксимації плоскої 
контактної поверхні набором трапецієвидних призм з вертикальними бічними сторонами (основою 1 вершиною). Сумарний 
гравітаційний вплив трапецій, що апроксимують контактні поверхні, з регіональним фоном, заданим ступеневим поліномом, 
визначається співвідношенням: 
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Для вирішення оберненої контактної задачі для поля, заданого похідною логарифмічного потенціалу на короткому ін- 
тервалі автором запропоновано систему лінійних інтегральних рівнянь зі швидкоспадними ядрами, дискретний вираз якої 
для загального випадку отримав назву узагальненого ітераційного процесу Лаврентьєва-Чорного 41: 
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Для аналітичного продовження гравіполя в верхній і нижній півпростір використовується методика Паркера на основі 
перетворень Фур'є згладженого вхідного сигналу. Для згладжування даних використовуються стандартні сплайнові методи і 
апроксимація методом спуску та найменших квадратів. Програмне забезпечення для перетворення картографічних коорди- 
нат перебуває в стадії розробки. Для вирішення систем лінійних рівнянь з розрідженими матрицями пропонується задіяти 
оптимальний метод біортогоналізації (5). Так, здійснюючи для системи лінійних рівнянь, до яких зводяться задачі геофізики 
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матриці А: К -о ро забезпечують нормальний перебіг процесу обчислень за рекурентними формулами 
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то розв'язок системи лінійних рівнянь Ах -Ф, увигляді Х - У. зі о; , отримаємо за п послідовних кроків наближень: 
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Ці програмні і математичні засоби органічно доповнять існуючу базу автоматизованих засобів інтерпретації даних 
потенціальних полів для вирішення задач як рудної так і структурної гравіметрії. 
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